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1 Abstract. 

' The data for the production of the ip' meson in pA collisions at 450 GeV 

^ ■ at CERN-SPS (of the NA50-collaboration) |I] yields a{il>' N) 8 mb under the 

assumption that the ij}' is produced as a result of the space-time evolution of a 
point-like cc pair which expands with time to the full size of the charmonium state. 
However, much higher values of (7{ip' N) are not ruled out by the data. We show 
' that recent CERN data confirm the suggestion of ref. 121 that color fluctuations 

' are the major source of suppression of the J/ip yield as observed at CERN in both 

pA and AA collisions. 
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1. Introduction 



, The NA50-collaboration observed quite different J/i/; and -^'-nucleon cross sections. 

' Analogous results have been found at Fermilab Both experiments found that 

, the charmonium-nucleon cross section is larger in the target fragmention region than 

at midrapidity. This is in good agreement with models that assume that the ■0' is 
produced as colorless, point-like cc pair, which expands with time to its full size and 
that there exists a relationship between the spatial distribution of color in a hadron 
O , and the cross section of its interaction with a nucleon. Such a relation is proved in 

pQCD (perturbative QCD) |^. For the nonperturbative regime it is a well known 
experimental fact that spatially larger hadrons have larger interaction cross sections. 
For example, 3.5 mb for the inelastic JZ-i/j-nucleon cross section was found at SLAC |E] 
while the inelastic 7r-nucleon cross section at this energies is 20 mb. From charmonium 
^ • models, e.g. in ref. |H| E] it is known that the different charmonium states {J /ip, X 

5^ , and if)') have different spatial sizes. We concluded in ref. [21 that these cross sections 

are dominated by non-perturbative contributions. 

The agreement of this scenario with the negative xp Fermilab data JH^ with was 
first demonstrated in ref. 0. Another attempt to describe these data with such an 
expansion scenario was done in ref. jj. However, in ref. @] 20%-50% of the suppression 
of charmonium states is due to a color octet state added to the model to explain also 
the large xp regime. This state as defined in ref. Q is in variance with QCD, because 
of the following reason. The eigen life time of this state adjusted to the Fermilab data 
is only 0.06 fm. A gluon emitted in such a short time has to have a momentum of 
1/(0. 06fm) K, 3.3 GeV relative to the cc pair j5|. In ref. ^ this gluon was assumed to 
be massless. 

In ref. it was shown that the production of J/ip^s in pA collisions can be 
described, if one takes into account the production and the subsequent decays of higher 
resonances (x,'0') into J/^'s. This leads to a significant increase of the absorption of 
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J/4>^s as compared to the propagation of genuine J/^-states. In ref. JUl it was shown 
also that the production of J/t/i's in AA cohisions is additionally suppressed by the 
final state interaction of charmonium states with newly produced particles like tt's, 
yo's and so on. 

In sect. 12 the used models are introduced and in sect. 13 they are compared to the 
data. A more detailed description can be found in ref. We summarize in sect.0] 



2. Model Description 

2.1. Semiclassical Glauber Approximation 

The suppression factor S for minimum bias pA collisions can be evaluated within the 
semiclassical approximation (cf. ^^1) as 

a{pA ^ X) 



a{pN X) 



J J d^B dz p{B, z) ■ exp ^- 



a{XN)p{B,z')dz'] . (1) 



A 

Here p{B,z) is the local nuclear ground state density (we used the standard 
parametrisation from (^{XN) is the interaction cross section of the charmonium 

state X with a nucleon and a{pA — > X) {a{pp X)) is the production cross section 
of the state X in a pA (pp) collision. We want to draw attention to the fact that this 
cross section changes with time due to the space-time evolution of color fluctuations. 
Therefore, it is necessary to keep a under the integral. The suppression factor S of 
J/i/j's produced in the nuclear medium is calculated as: 

= 0.6 • (0.92 • S''/''' + 0.08 • S"^'') + 0.4 • . (2) 

Here S'^ are the respective suppression factors of the different pure charmonium states 
X in nuclear matter. Eq. ||2Jl accounts for the decay of higher resonances after they 
left the target nucleus into J/4>^s. The fractions of J/4>'s that are produced in the 
decays of higher resonances in eq. are taken from ref. However, in ref. |2I it 
is assumed that the different charmonium states interact with nucleons with the same 
cross section, which is in disagreement with the data from the refs. 0E|- 

In line with the above discussion we want to stress here that Eqs. ^ and (|2Jl 
are applicable at CERN energies for central and negative rapidities, but have to be 
modified, if applied already at r/cm. ^ at RHIC or higher energies, because at 
higher energies charmonium states can be produced outside of the nucleus and the cc 
pairs propagate through the whole nucleus without forming a hadron. Data are often 
presented in the form a vs. A, where a is given by 

apA = (Tpp ■ A" so that 5* = A"-^ . (3) 

In the semiclassical Glauber approximation, we take into account these color 
fluctuations in an effective way as described in ref. . We assume that charmonium 
states are produced at z as small cc configurations predominantly through gluon-gluon- 
fusion, then they evolve - during the formation time - to their full size. Please note 
that there is up to now no theoretical or experimental proof for the assumption that 
charmonium states are produced in point-like configurations as predicted in pQCD. A 
way to test this experimentally was suggested recently in ref. It's based on the 
idea that the -0' and the ip" are mixed states of the 2S and the ID charmonium states. 
The names of these states comes from a comparison of nonrelativistic charmonium 
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models with the nonrelativistic wave functions of the positronium (actuaUy even the 
name charmonmm was given to pronounce the similarity with positr omwrn). Namely 

IV'') = cose* |2S') + sin 61 III?) , 

IV-") = cos6l 1 1Z3) - sin 61 1 25') . (4) 

Since only the S'-wave contributes to the decay of ip states into e+e^-pairs (at least 
in nonrelativistic charmonium models) the value of 6* — 19 ± 2° can be determined 
from the data on the e+e" decay widths of ip' and -0". If the production of these 
charmonium states is pointlike, then only the S'-wave is produced. These leads to an 
universal ratio of 

a(V'') r(V'' ^ i+i-) ■ 

<j{ip") resp. (j{ip') are here the production cross sections of the tp" and resp. tp' in 
various processes. Predictions for different processes can be found in ref. 

If the formation length of the charmonium states, If, becomes larger than the 
average internucleon distance (If > r^jsf « 1.8 fm), one has to take into account the 
evolution of the cross sections with the distance from the production point ^S] . Here 
we assume motivated by quantum diffusion that the cross section increases linearly 
with time. The formation length of the J/ijj is given by the energy denominator 
If « i i — , where p is the momentum of the J/ip in the rest frame of the target. 

With p ~ 30 GeV, the momentum oi a, J/ip produced at midrapidity at SPS energies 
(Eiab ~ 200 AGeV), this yields // « 3 fm, i.e. a proper formation time of r/ = 0.3 
fm. As formation time of the tp' in its rest system we use here the radius given by 
nonrelativistic charmonium models, e.g. see the refs. 0111. This radius is r = 0.45 
fm for the ip'. A larger value of t/ for the ip' is supported also by the extraction of 
the formation time of the J/ip |17j . Finally the formation time is tf = J ■ Tf, where 
7 is the the Lorentz- factor of the charmonium state relativ to the nuclear target. For 
higher gamma factors, i.e. at higher energies, the formation time becomes larger than 
the nuclear targets. In this regime a hadronic description during the formation time 
is questionable. A partonic model for this energy range was proposed in ref. |18|. 



2.2. Vector Dominance Model for jA 

The VDM (Vector Dominance Model) takes into account only the direct diffractive 
production of the J/tp and the ip' , while the GVDM (Generalized Vector Dominance 
Model) accounts also for the non-diagonal transitions {ip' + N J/ip + N and 
J/ip + N — > Ip' + N). The later are needed, because in photoproduction the 
particles are produced as point like configurations and develop then to their average 
size. In a hadronic model like the GVDM this is taken into account in form 
of the interference due to the nondiagonal matrix elements. In the GVDM the 
photoproduction amplitudes f~f^ and f-^^i for the J/ip and the ip' are given by |19| 

/71/1 — f^4> H 7 /«/''»/' 

6 6 

f-ii,' = -rf^i>' + -^f^p'Tp' (6) 
U U' 

Here and are the J/ip ~ 7 and the ip' — ^ coupling and fvv are the amplitudes 
for the processes V+N — * V' + N, where V and V' are the J/ip and the ip' respectively. 
In the VDM the non-diagonal amplitudes with V V are neglected. The importance 
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of the nondiagonal transitions is evident, because the left hand side of eq. ^ is small. 
If it is neglected as a first approximation |19j . then /^'^ = —j^fiptp ~ 1-7 • /^^. And 
due to the CPT-theorem /^/^ = /^^'. 

In ref. [20] the GVDM yields approximately 8±2 mb for the i/j'-nucleon interaction 
cross section at SPS-energies as can be seen in Fig. 1. is the laboratory energy of 
the photon. The J/i/i-nucleon interaction cross section at SPS-energies approximately 
3.5 - 4 mb is used as input into the analysis of ref. UHl- The accuracy of such GVDM 
in predicting the t/i'N cross sections is not clear. The above calculation demonstrates 
that implementing color transparency leads to significantly larger cross sections of the 
tp' N interaction. 
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Figure 1: The energy dependence 
of the elementary charmonium- 
nucleon cross sections found in the 
GVDM. The filled areas show the 
variation of the cross sections due to 
the uncertainty of the experimental 
J/t/jN cross section. 



3. Comparison with NA50 and NA51 Data 

In Fig. 2 we show a comparison between calculations with different cross sections 
and different expansion times and the NA50 data ^ for pA collisions and the NA51 
data [2] for pp and pD collisions for the cross section of V'' N interaction vs. the mass 
of the target. The y-axis shows B^^^a^i /A where -B^^ is the branching ratio for the 
decay of the ip' into dimuon pairs, and a^i is the production cross section. The "5.1 
mb, instant formation" curve in Fig. 2 is the fit of the NA50 collaboration to their 
data. Instant formation means that they assumed that the ip' is produced with the 
full cross section and not as a point like particle as in the description of this paper. 
(Note the NA50 collaboration fitted B^^a^,i /(Jdy, where DY means Drell-Yan, we 
multiplied this fit with the DY cross section in pp collisions measured by NA51). 

The "8 mb, i/ = 0.45 fm" curve is the eye-ball fit of the model described in this 
paper. For the comparison with the data we need the production cross section of the 
■!/;' in pp collisions as input. We used here the average of the pp and pD data of the 
NA51 collaboration. The value of 8 mb agrees well with the model parametrisations 
discussed in the sections |2J However, we compare also with the calculation with the 
parameters of ref. 0, i.e. a{ip'N) — 20 mb and tf ~ 0.6 fm. For this comparison 
we used the production cross section of the tp' in pD collisions divided by two as 
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input. This is also close to the value of the NA50 fit. One can see in Fig. 2 that the 
calculation with these parameters is also in good agreement with the data. A value for 
a{'4''N) of the size of 20 mb is favored by the nucleus-nucleus data as shown in ref. [TU] 
and in Fig. 3. Plotted is Bfj_fj_a{J/ip)/a{DY) with the absorption cross sections from 
ref. and the NA50 data 122 for PbPb collisions vs. the transverse energy Et, a 
measure for the centrality of the collision. The calculation agrees well with the data. 
The calculation for the V'' underestimates the data. However it is not understood, if 
this is due to the high value of a{ip'N) — 20 mb, or if nondiagonal transitions like in 
sect. 12.21 should be taken into 
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Figure 2: Plotted are the -Bp^cr^/ /A 
values extracted from calculations 
with different absorption cross sec- 
tions and different formation times, 
the NA51 data for pp and pD, and 
the NA50 data for pBe, pCu, pAg, 
and pW vs. the mass A of the tar- 
get. 
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account in A A collisions, too. 

The value of 8 mb is smaller than the theoretical estimate 20 mb of ref. 2 . This 
is because in ref. (2j a formation time of 0.6 fm was chosen for the ^' , while we used 
her 0.45 fm, the radius of the tp' given by the charmonium models. The fact that 
the formation time is not known very well is another theoretical uncertainty. Further 
uncertainty comes from using diffusion model of expansion at the distances comparable 
to the scale of the soft interaction. Within the error bars the ^'-nucleon cross section 
extracted from these pA data and the predictioir of the GVDM, discussed in sect. 12.21 
are qualitatively similar. However, further data are needed to learn more about this 
cross section. 
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Figure 3: Plotted is 

Bf,f,a{J/ii)/cj{DY) with tlie 
absorption cross sections from 
ref. 12 and the NA50 data for 
PbPb collisions vs. the transverse 
energy Et 



4. Conclusions 

The new data of the NA50-collaboration Tl and the data of the E866-collaboration [3] 
prove that the ip'-mxcleon cross section is much larger than the J/^-nucleon cross 
section. This is in agreement with the photoproduction data for these charmonium 
states as discussed in the framework of the GVDM in section 12.21 This confirms 
the QCD prediction that the strength of hadron-hadron interactions depends on the 
volume occupied by color. 

Within the assumption that charmonium states are produced as point like white 
states, we demonstrated that the data jlj can be fitted with a V^'-nucleon cross section 
of a{ip'N) w 8 mb. However, a much larger cross section of e.g. a{ip'N) « 20 mb is 
not ruled out by the data. Due to the large experimental errors we conclude that the 
data and the QCD-motivated models agree, but further data with higher accuracy 
and covering larger rapidity range are needed. 

This cross section will be measured soon in proton-nucleus collisions at HERA B 
at an energy of Eiab = 920 GeV. The advantage of this experiment is that it covers 
a larger range of Feynman-XF, especially in the negative xp region. In this region 
effects due to the formation time of the hadron will be less important and the genuine 
cross sections will be measured. 
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